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X-ray investigation of the diffuse emission around plausible 7-ray 
"p^". emitting pulsar wind nebulae in Kookaburra region 

O ■ Tetsuichi Kishishita^, Aya Baniba^, Yasunobu Uchiyama^''*, Yasuyuki Tanaka^, and 

. ■ Tadayuki Takahashi^'^ 

ABSTRACT 

(N 

We report on the results from Suzaku X-ray observations of the radio com- 
W . plex region called Kookaburra, which includes two adjacent TeV 7-ray sources 

HESS J1418-609 and HESS J1420-607. The Suzaku observation revealed X-ray 
diffuse emission around a middle-aged pulsar PSR J1420-6048 and a plausible 
PWN Rabbit with elongated sizes of ax = 1'.66 and ax = 1'.49, respectively. 
The peaks of the diffuse X-ray emission are located within the 7-ray excess maps 
j^ . obtained by H.E.S.S. and the offsets from the 7-ray peaks are 2'. 8 for PSR J1420- 

6048 and 4'. 5 for Rabbit. The X-ray spectra of the two sources were well repro- 
duced by absorbed power-law models with F = 1.7 — 2.3. The spectral shapes 
Q\ . tend to become softer according to the distance from the X-ray peaks. Assum- 

^ . ing the one zone electron emission model as the first order approximation, the 

in '• ambient magnetic field strengths of HESS J1420-607 and HESS J1418-609 can 

cn . be estimated as 3 /iG, and 2.5 /iG, respectively. The X-ray spectral and spa- 

^ . tial properties strongly support that both TeV sources are pulsar wind nebulae, 

in which electrons and positrons accelerated at termination shocks of the pulsar 
winds are losing their energies via the synchrotron radiation and inverse Compton 
/\ '. scattering as they are transported outward. 
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Subject headings: gamma-rays: observations — ISM: individual (Rabbit) 
stars: pulsars: individual (PSR J1420-6048) — X-rays: ISM 



INTRODUCTION 



Gamma-ray emitting Pulsar Wind Nebulae (PWNe) are excellent test grounds for study- 
ing pulsars' relativistic winds and particle acceleration that takes place presumably at termi- 
nation shocks. Especially, PWNe which accompany the diffuse synchrotron X-ray emission 
show the important clues to understand the time evolution of the PWNe and how the ac- 
celerated particles escape from the shocks. While young PWNe with characteristic ages 
of less than 10 kyr are relatively well studied in the past obs ervations with Chandra and 
XMM-Newton (JKargaltsev fc Pavlovl l2008l : iBamba et al.ll2010al ). the sample of middle aged 
PWNe remains small. The Kookaburra region includes a middle aged PWN (PSR J1420- 
6048) and a plausible PWN (Rabbit), making it a suitable target for the detailed study of 
the diffuse emission with the Suzaku X-ray observatory. Suzaku, characterized by the low 
detector background compared to Chandra and XMM-Newton, is crucial for the analysis of 
faint and diffuse X-ray emission. 

The complex of compact and extended radio/X-ray sources, called Kookaburra (desig- 
nated by [Roberts et al.lll999l ). spans over about one square degree along the Galactic plane 
around / = 313°. 4. The H.E.S.S. galactic survey revealed two very high energy (VHE) sources 
in this region; the brighter of the two sources, HESS J1420-607, is centered at the position 
of (R.A., Dec.) = (14^20'"09^ -60°45'36") with an intrinsic extension of axev = 3'.3 ± 0'.5, 
and the slightly less bright second source, HESS J1418-609, is centered at the position of 
(R.A., Dec.) = (14^18'^04^ -60°58'31") with an intrinsic extension of cTTev = 3'.4 ± 0'.6 (a 
major-axis of 4'. 9 ± 1'.5 and a minor-axis of 2'. 7 ± 0'.7 fitted with an elongated Gaussian 
shape). The Kookaburra region is a lso bright in the GeV band and registered as a Fermi 
7-ray source (IWeltevrede et al.ll2010[ ). 



There is an energetic pulsar PSR J1420-6048 at the south of HESS J1420-607. It is 
located at (R.A., Dec.) = (14^^20™08^20, -60°48'17".2), which is ~ 3M offset from the 
center of the VHE 7-ray emission. PSR J1420-6048 is a radio /X-ray pulsar with a period o f 
P = 68 ms and period derivative of P = 8.3 x 10"^^ jD'Amico et aUboOsI : iNg et al.ll2005[ l. 
The distance of the pulsar was estimated as d = 5.6 kpc ba sed on the pulsar's disp ersion 



measure with the NE2001 Galactic electron-density model ( ICordes &: Lazio 1 12002[ ). The 
characteristic age and the spin-down luminosity are r = 13 kyr and E = 1.0 x 10^^ erg s~^, 
respectively. The strength of surface magnetic field is Bs = 2.5 x 10^^ G. 
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Chandra obse rvations revealed a faint and narrow diffuse X-ray emission in K3 neb- 
ula (designated by [Roberts et al.lll999l with the radio observation) around PSR J1420-6048 
( JNg et al.ll2005l ). Spectral fitting with an absorbed power-law model to the K3 nebula (2' 



radius from the pulsar) yielded an absorption column density of A^h 
with a photon index of F 
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The Rabbit nebula is located at the Eastern edge o f HESS J1418-60 9 with a distance 
of 8'. 2 to the best central position of the VHE emission. iNg et al.l ( 120051 ) found two point- 
like sources with the XMM-Newton data, labeled as Rl and R2; the brighter source (Rl) 
is located at the edge of the Rabbit, while the fainter source (R2) at (R.A., Dec.) = 
(14'^18™39*'.90, -60°57'56".5), which appears embedded in the narrow diffuse emission. Tim- 
ing analysis with the EPIC pn data revealed a P = 108 ms pulsation with a period derivative 
of P = 1.07 X 10^^^ from R2, although not highly significant. Assuming these values are 
correct, the characteristic age can be estimated as r = 1.6 kyr. Fitting the spectrum from 
the whole Rabbit nebula (3' radius) with an absorbed power-law model, they found an ab- 
sorption column density of A^h = 1.4 ±0.2 x 10^^ cm~^ and a photon index of F = 1.5 ±0.14. 



The two VHE sources are considered to be PWNe (JRoberts et al.lll999l : iNg et al.ll2005l ). 
given the X-ray and TeV results reported so far. Generally, detailed observations of the 
synchrotron X-ray nebulae are crucial to understand spatial and spectral distributions of 
accelerated electrons. However, previous X-ray observations did not have enough sensitivity 
to detect faint X-ray nebulae with sizes comparable to the TeV nebulae. In this paper, we 
report new detections of the diffuse X-ray emission of the Kookaburra complex using the 
Suzaku satellite. 



OBSERVATION AND DATA REDUCTION 



We observed the Kookaburra region with the Suzaku satellite in 2009 F ebruary. The 
Suza ku observations were performed with the X-ray Imaging S pectrometer (XIS: iKovama et al, 



20071 ) in 0.3 - 12 keV and the Hard X-ray Detector (HXD: iTakahashi et al 



600 k eV. The XIS, located at the focal plane of the X-ray telescopes (XRT: ISerlemitsos et al. 



m 



20071) in 13 - 



20071 ). consists of one back-illuminated CCD camera (XISl) and three front-illuminated 
CCDs (XISO, 2, and 3). One of the front-illuminated CCDs, XIS2, was not available at the 
time of our observations, since it suffered from a fatal damage on 9 November 2006, and 
unusable since then. The XIS instruments were operated in a nor mal full-frame clockiii g 
mode (a frame time of 8 s) with Spaced-row Charge Injection (SCI) ( iNakajima et al.ll2008l ). 
The HXD consists of the silicon PIN photo diodes (hereafter PIN) capable of observations 
in the 13 - 70 keV band and the GSO crystal scintillators (hereafter GSO) which cover the 
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40 - 600 keV band. Since we could not find significant signals in the HXD data after sub- 
tracting the non- X-ray background (NXB), cosmic X-ray background (CXB), and galactic 
ridge components, we focused on the XIS data analysis in this paper. 

We used data sets processed by a set of software of the Suzaku data processing pipeline 
(version 2.2.11.24). The telemetry saturating time was excluded in the pipeline processing. 
Basic analysis was done using the HEASOFT software package (version 6.4). We made use 
of cleaned event files, in which standard screening was applied. The standard screening 
procedures include event grade selections, and removal of time periods such as spacecraft 
passage of the South Atlantic Anomaly (SAA), intervals of small geomagnetic cutoff rigidity 
(COR), and those of a low elevation angle. Specifically, for the XIS, the elevation angle larger 
than 5° above the Earth and larger than 20° from the sunlit Earth limb are selected. The 
total exposure from the two observations amounts to 72 ks after standard data screening. 
Table [1] gives the log of the Suzaku observations. 

Table 1: Summary of the Suzaku observation on the Kookaburra region 



Obs. ID 



Coord. (J2000) 
R.A., Dec. 



Exposure 
XIS 



Date 



503110010 215.0292,-60.8156 50.3 01/11/2009 
503071010 214.6625,-60.9675 21.3 02/14/2009 



° Units are ks. 



3. ANALYSIS AND RESULTS 



3.1. Image Analysis 

In order to combine two different pointing images, we extracted the photon events in 
the energy range of 2-10 keV from each sensor. The data between 5.73 and 6.67 keV were 
removed from the image to exclude the calibration sources. We corrected the vignetting 
effect by dividing the image wit h a flat sky image si mulated in the energy range of 2-8 keV 
using the XRT + XIS simulator (llshisaki et al.ll2007l ). The image was binned to 8 x 8 pixels 
and smoothed with a Gaussian function of 3(j. A combined Suzaku XIS (0+1+3) image in 
the Kookaburra region is show in Fig. [H Both HESS J1420-607 and HESS J1418-609 
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are reported as diffuse TeV-7-ray emitting sources (jAharonian et al.ll2006l ). The 
intrinsic extensions of the diffuse emissions are indicated as white circles in the 
figure. The overlaid contours are taken from the H.E.S.S. observation. We can 
clearly see that the separate X-ray sources locate with some offsets from the 
7-ray peak positions and are coincident with the VHE sources. 



K3/PSR JI42O-6O48 



We extracted the photon events in the energy ranges of 1-2 keV and 2-10 keV. The XIS 
images are shown in Fig. [2J We can see several thermal sources in the lower energy image, 
which are invisible in the higher energies. The bright X-ray emission is well coincident with 
the location of PSR J1420-6048 and extended diffuse emission surrounds the bright pulsar. 
The diffuse component has an elongated shape which extends from the pulsar to the 7-ray 
peak. In order to determine the extension of non-thermal X-ray emission, we 
created a surface brightness profile from the enclosed region shown as "rect" 
in Fig. [2] (a) along the north to south direction. The surface brightness in 
2—10 keV vs. relative coordinate is shown in Fig. [Si The profile was fitted 
with a Gaussian function plus constant as a background to derive the size of the 
diffuse emission, using the sigma of the Gaussian (fix)- Since a Gaussian profile 
often provides a rea sonable approximation for the surface brightness profiles 
( jBamba et al.l l2010b ). we used relative coordinates between 2'. 2 and 9'.0 as a 
fitting range in consideration for the non-symmetrical profile. We excluded a 
range between 2'. 5 and 4'.0 as a narrow pulsar component. This is to avoid the 
contribution from the bright pulsar, whose image is smeared by the point-spread 
function of XRT. We defined the source size as three times the a^- The fitted 
result shows crx = 1'.66 ± 0'.09. If we assume the distance as 5.6 kpc, the physical 
size of the extended diffuse emission corresponds to 8.1 ±0.4 pc. The intrinsic 
emission size of the bright pulsar was also estimated by fitting the image with 
a two-dimensional Gaussian function, however, the calculated size wa s smaller 
than the attitude fluctuation of the satellite (see lUchiyama et al.ll2008l ). 



Rabbit 



Fig. a shows the XIS images of the Rabbit nebula. Similar to K3/PSR J1420-6048, 
thermal sources are visible in the lowe r ener gies. The X-ray peak is coincident with the 
point-like sources reported by iNg et al.l ( 120051 ) . which are indicated as black cross points in 
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Fig. 1. — Suzaku XIS (0+1+3) image around the Kookaburra region in a linear scale. The 
extracted energy range is 2- 10 keV. Green contour s denote the 5 a, 7.5 a and 10 a significance 
levels reported by H.E.S.S. ( lAharonian et al.ll2006h . The white indicates the best-fit positions 
and intrinsic extensions of the TeV-7-rays. 
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Fig. 2. — XIS images around K3/PSR J1420-6048 in the energy range of (a) 1-2 keV and (b) 
2-10 keV. PSR J1420-6048 is depicted a s a black diamond. The white cross point corresponds 
to the best-fit position of VHE 7-rays ( lAharonian et al.ll2006l ). 
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Fig. 3. — Projected profile of the surface brightness around K3/PSR J1420-6048. The 
extracted region is shown as the white rectangle ("rect") in Fig. |2](a). 



Fig. m According to iNg et al.l (|2005[ ). these bright sources are embedded in narrow diffuse 
emission. Compared with the XMM-Newton, the diffuse emission in Fig. H] (b) looks more 
extended than the previous image, and another diffuse structure which is clearly separated 
from the point-like sources can be seen in the southeast direction. In order to determine 
the size of the diffuse emission, we made a projected profile with the same procedure as 
in K3/PSR J1420-6048. The extracted area is indicated as "rect" in Fig. 1 (a). The ID 
projected result is shown in Fig. [51 We excluded a range between — 1'.2 and — 0'.4 from 
the fitting as the central pulsar region to avoid its influence. The projected profile was well 
represented with a Gaussian function plus constant background component. The size of the 
diffuse component is estimated as ax = 1'.49 ± 0'.07 from the ID profile. The physical size 
of the diffuse emission can be estimated as 6.5 ± 0.3 pc. We assumed the distance as 5 kpc. 



3.2. Spectral Analysis 

K3/PSR JI42O-6O48 

Firstly, we checked the background level extracted from the "bgd" region in Fig. [2] (b). 
We compared its background spectrum with public blank sky observations on Lockman-hole 
(Obs. ID: 102018010) and also with ASO 0304 (Obs. ID: 504054010), which locates near 
the Kookaburra region, i.e. (R.A., Dec. ) = (213. 3355, -62.0808). The background levels are 
consistent within 10% and the Galactic ridge emission is ignorable. We thus concluded that 
the "bgd" region is a source-free in the following spectral analysis. 

In order to determine the absorption column density around PSR J1420-6048, we ex- 
tracted the photon events within an angular distance of 1'.5 from the pulsar, which corre- 
sponds to "ringl+ring2" in Fig. |2] (b). We have co-added the data from XISO and XIS3 
to increase statistics. The response (rmf) files and the auxiliary response (arf) files were 
produced using xisrmfgen and xissimarfgen, respectively. Fig. E] shows the spectral result 
with the joint-fitted XISO+3 (black) and XISl (red) data. The spectral shapes are well 
represented by an absorbed power-law model ('wabsxpow'). From the spectral fitting, the 
absorption column density was estimated as Nh = 4.4 ±0.3 x 10^^ cm~^, which is consistent 
with the previous observations. To determine the total flux from the pulsar and diffuse 
emission, we used the "regl-reg2" as the extracted region. The flux in the energy range of 
2-10 keV is 4.56 ± 0.10 x 10~^^ erg-cm~^ • s~^. The systematic uncertainty coming from the 
difference of the effective area between the arf responses for diffuse and point sources was 
less than 10%. 

In order to investigate the spatial dependence of the spectral shapes, we separated the 
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Fig. 4. — XIS images around the Rabbit nebula in the energy range of (a) 1-2 keV and 
(b) 2-10 keV. The black cross points show the point-like sources detected by XMM-Newton, 
denoted as "Rl" and "R2" in iNg et al.l (120051) . The w hite cross point corresponds to the 
best-fit position of VHE 7-rays (jAharonian et al.ll2006l ). 
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Fig. 5. — Projected profile of the surface brightness around the Rabbit nebula. The extracted 
region is shown as the white rectangle ("rect") in Fig. H] (a). 
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extraction region into four annular regions, i.e. "ringl": 0-0'. 5, "ring2": 0'.5-l'.5, "ringS": 
l'.5-2'.5, "ring4": 2'. 5-3'. 5. To reduce the uncertainty of the column density with limited 
statistics, we fixed the absorption column density at Nh = 4.4 x 10^^ cm~^ determined from 
the "ringl+ring2" region. Fig. [7] shows the spectra (XISO+3) for each circular region. In 
order to exclude the thermal contamination, we used the 2-10 keV band for the spectral 
fitting. In the "ringS" spectrum, we can see a weak trend of the iron emission line around 
6.7 keV. This might come from unresolved thermal sources. The photon index vs. relative 
coordinate is shown in Fig. |SJ Quoted errors are at the la confidence level. We can see that 
the photon index slightly changes according to the distance from the bright pulsar region. 
The contamination effect from a point source in "ringl" on the outer region "ring4" was 
estimated as ~5%. The fitted parameters are summarized in Table [2l As for the HXD-PIN, 
we could not detect any significant pulse in timing analysis, thus there is no excess in the 
hard X-ray band from the pulsar. 



Rabbit 

We performed the same analytical procedures to the photon events from the Rabbit 
nebula as described in K3/PSR J1420-6048. Fig. |9] shows the spectra extracted within an 
angular distance of 1'.66 from the X-ray peak, which corresponds to the "ringl+ring2" region 
in Fig. m (b). The background spectrum was consistent with that of the blank sky data. 
The spectral shapes are well represented by an absorbed power-law model. The absorption 
column density was determined as Nh = 2.7 ± 0.2 x 10^^ cm~^ by the spectral fitting. 

For the spatial dependence of the spectral shapes, we chose the annular regions as 
"ringl": 0-0'.83, "ring2": 0'.83-l'.66, "ring3": l'.66-3'.0, "ring4": 3'.0-5'.0. The spectral 
results are shown in Fig. [TOl To exclude the thermal contamination, we used the energy 
range of 2-10 keV for spectral fittings. The photon index vs. relative coordinate is shown in 
Fig. [TTl The spectral shape becomes softer according to the distance from the inner bright 
core to the outer region, which is a similar trend as in the K3//PSR J1420-6048 region. The 
systematic error of the PSF in "ringl" to "ring4" was less than 5%. The best-fit parameters 
are summarized in Table |2l The total flux in the energy range of 2-10 keV is estimated as 
6.27 ± 0.13 X 10^^^ erg-cm~^ • s~^ determined from the extraction region of "regl-reg2". 
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Fig. 6.— Spectra of the "rmgl+ring2" in K3/PSR J1420-6048 (black: XISO+3, red: XISl). 
Solid lines show the fitted absorbed power-law model. The bottom panels show residuals 
from the best-fit. 
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Fig. 7. — Spectra of the circular regions in K3/PSR J1420-6048. The extracted regions are 
shown in Fig. |2] (b). 
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Fig. 8. — Spatial dependence of the photon index from the bright pulsar PSR J1420-6048. 
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Fig. 9.— Spectrum of the "rmgl+ring2" in Rabbit (black: XISO+3, red: XISl). Solid 
lines show the fitted absorbed power-law model. The bottom panels show residuals from the 
best-fit. 
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Fig. 10. — Spectra of the circular regions in the Rabbit nebula. The extracted regions are 
shown in Fig. H] (b). 
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Fig. 11. — Spatial dependence of the photon index from the bright center region in the 
Rabbit nebula. 
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Table 2: Results of Suzaku XIS Spectral Fitting 



Object & region 


N^ (wabs) ^• 


^Photon index "^ 


Flux [2-10 keV] ^' 


' X» 


K3/PSR J 1420-6048, ringl+ring2 


4.4 ± 0.3 


2.00 ± 0.10 


1.53 ± 0.03 


1.01 (72) 


regl-reg2 


4.4 


(fix) 


2.17 ±0.06 


4.56 ±0.1 


0.83 (33) 


ringl 




— 


1.69 ±0.16 




0.86 (32) 


ring2 






2.06 ±0.06 




1.01 (32) 


ringS 






2.28 ±0.07 




0.86 (35) 


ring4 




— 


2.29 ±0.09 




0.94 (35) 


Rabbit, ringl+ring2 


2.7 ± 0.2 


1.82 ± 0.10 


2.65 ± 0.05 


0.87 (101) 


regl-reg2 


2.7 


(fix) 


2.00 ± 0.06 


6.27± 0.13 


0.80 (52) 


ringl 






1.77 ±0.08 




1.06 (46) 


ring2 




— 


1.94 ±0.06 




1.09 (35) 


ringS 




— 


2.12 ±0.06 




1.10 (46) 


ring4 






2.21 ±0.08 




1.08 (33) 



"Quoted errors are at the \a confidence level. 

^Nh units are 10^^ cm^^. 

^Unabsorbcd fluxes are given in units of 10^^^ erg cm^^ s^-'^. 
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4. DISCUSSION 

4.1. Origin of the X-ray and 7-ray Emission 

We have detected the diffuse X-ray emission around K3/PSR J1420-6048 and Rabbit, 
which were unnoticed in the previous observations, in the Kookaburra region with the Suzaku 
sateUite. The X-ray peaks of the two sources are both within the error circles of the TeV 
sources and each of them has an offset from the 7-ray peak with 2'.8lo;g ^r K3/PSR J1420- 
6048 and 4'.5l;Q;g for Rabbit. The XIS spectra of the two sources were well reproduced by 
an absorbed power-law model with a photon index of F = 1.7 - 2.3. 

If TeV/X emission comes from the same object, the origin of the VHE 7-rays can 
be explained by the inverse Compton (IC) scattering of the CMB photons by high-energy 
electrons, while the X-ray emission via synchrotron radiation in a mean magnetic field B. In 
this case, typical energies of responsible electrons are i^g — 20 ■ {y^j) TeV for IC 7-rays 
(in the Thomson regime) at a photon energy eic, and i^e — 70 • (yqIIg)" • (rkeV) TeV 
for X-rays at esyn, respectively. The synchrotron and IC spectra produced by the relativistic 
electrons obeying an energy distribution of Ne{E) oc E~^, have the same spectral shape, 

-^sync or icl^syn or ic) OC Cgy^^^j. j^. as loug as IC Scattering occurs in the Thomson regime. If 
the high-energy electrons are continuously injected into the radiation zone at a constant rate 
and loosing energy predominantly by synchrotron or IC (Thomson) losses, the electron energy 
spectrum becomes one power of E steeper above the cooling break E]^,. that is determined by 
the cooling timescale and the age of the source. As a result, the synchrotron and IC spectra 
from the cooled electrons become softer by AF = 0.5. 

Our Suzaku results clearly show that the photon index smoothly changes according 
to the distance from the bright center regions, and the X-ray photon indices near the 7- 
ray peaks are roughly consistent with those in the TeV energies in both objects (Ftbv 



2T7± 0.1 for HESS J1420-607 and Fxev ^ 2.22 ±0.1 for HESS J1418-609, lAharonian et al. 



( I2OO6I )). In addition to that, the difference of the photon index between the bright center 
and diffuse emission is AF ~ 0.5. The smooth spectral steepening of the synchrotron X- 
ray emission likely reflects synchrotron burn- off of the accelerated elec trons. Such effects 



have been observed for some young P WNe (JGaensler and Sland |2006| ) , but we now find 



an interesting example of the smoothly steepening of the synchrotron X-ray emission in a 
middle-aged PWN of K3/PSR J1420-6048 (rsd ~ 13 kyr). 

In order to estimate the mean magnetic field strength, we modeled the spectral energy 
distribution. Since information about the spatial dependence of the TeV emission is not 
available, we assumed a simple one-zone synchrotron + IC model as the first order approxi- 
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mation, in which a single population of relativistic electrons in TeV energy range emit both 
X-rays through synchrotron radiation and TeV 7-rays through IC scattering off the CMB. 
The electrons energy distribution is formally assumed to be an exponentially cutoff power- 
law of the form N^. oc 7~Pexp(— 7/7m), where 7 = Ee/nieC^. We do not specify the physical 
meaning of 7m, which may account for both the maximum electron energy and the cooling 
break. This spectral form is invoked only for the purpose of an estimate of the magnetic 
field. 

Fig. [12] shows the spectral energy distribution of HESS J1420-607 and 
HESS J1418-609 at the X-ray and 7-ray bands; the model curves are obtained 
with a combination of 7^ = 8x10^ (HESS J1420-607) or 7^ = 2x10® (HESS J1418- 
609) and p = 2.3 for both objects. The black thick lines overlaid on the Suzaku 
data points indicate the absorbed power-law models with best-fit parameters as 
shown in Table [2l We note that the "ringl" region mainly comes from bright 
pulsars while the "regl-reg2" region includes both emissions from pulsar and 
diffuse components. The deviation between the spectrum of the "regl-reg2" 
region and the model curve in higher X-ray energies comes from the pulsar's 
intrinsic hard spectrum since the intrinsic diffuse emission can be approximately 
estimated as "regl-reg2" - "ringl". The magnetic fields can be estimated as 3 fiG for 
HESS J1420-607 and 2.5 /iG for HESS J1418-609. These values are in agreement with the 
magnetic field strength that is expected for middle-aged PWNe, which lends support to the 
PWN scenarios for both objects. 



4.2. Morphological properties 

The diffuse X-ray emission of K3/PSR J1420-6048 has asymmetry (a circular shape of 
the radius i?diffuse(3crx) = 8.1 ± 0.4 pc + a tail-like component) and its tail-like component 
extends to the 7-ray peak. On the other hand, the diffuse emission in the Rabbit seems to be 
more concentrated near the bright central source with higher surface brightness and relatively 
symmetrical shape of i?diflfusc(3crx) = 6.5 ±0.3 pc. The angular sizes of the VHE 7-ray sources 
are larger than those of the X-ray emission regions, e.g. -Rdiffuse(3crTev) = 16.1 ± 2.4 pc 
(HESS J1420-607) and /2diffusc(3aTev) = 14.8 ± 0.9 pc (HESS J1418-609). 

For the estimated magnetic field of S ~ 3fiG, the energy loss rate due to synchrotron 
radiation is just comparable to the loss rate via IC scattering off the CMB in the Thomson 
regime, which is given by 

--~140(^)"ky.. (1) 
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The lifetime of the electrons emitting synchrotron X-rays at 1 keV is Tsyn(e = IkeV) ~ 
10 kyr for the estimated magnetic field strength. On the other hand, the energy loss 
time scale of the electrons emitting IC 7-rays at 1 TeV is ric(e = 1 TeV) ~ 70 kyr (see 
e.g. 



de Jager &: Djannati-Ataill2008l ). The TeV emitting electrons have a longer lifetime 



than the cooling time of the X-ray emitting electrons, which is shorter than the characteris- 
tic age of PSR J1420— 6048 and probably also shorter than that of the pulsar in the Rabbit. 
Therefore the morphological differences between th e two bands can be expected, though 



details may depend largely on the PWN evolution ( jBamba et al.l l2010bl : IVink et al.l 12011 
Van Etten k Romanill201ll ). 



According to iMattana et al.l (|2009[ ). 7-ray to X-ray energy flux ratio of PWNe is pro- 



portional to the characteristic age. The energy flux ratios for both objects are -^^(l — 
30 TeV)/Fx(2-10 keV) ~ 3.2 for PSR J1420-6048 and F^(l-30 TeV)/Fy(2-10 keV) ~ 0.16 



for the Rabbit nebula (-^7(1 — 30 TeV) values were taken from iMattana et al.l (120091 )). If 



we assume that the characteristic ages for both objects are correct, i.e. r = 13 kyr for 
PSR J1420-6048, and r = 1.6 kyr for Rabbit, the flux ratios roughly follow the Mattana's 
relation. 



4.3. Diffusion of Electrons 

The high-energy electrons are thought to be transported via convection and/or diffusion 
due to the magnetic fields randomized at the termination shock. Let us assume that diffusion 
is the dominant transport mechanism. Indeed the pulsar wind slows down from relativistic to 
relatively low expansion velocities at some distance after the pulsar wind has been shocked. 
The diffusion coefficient can be estimated as 



K 



R: 



/Ts, 



■diffuse/ 'syn, 



(2) 



where -RdiSuse is the observed size of the synchrotron X-ray nebula. On the other hand, the 
diffusion coefficient can be written as 



Ac/3 

2.6 X 10^^ ■ / • 



B 



10 fiG 



E, 



70 TeV 



cm^ s ^ . 



(3) 
(4) 



where A = /rL is the mean free path of scattering, parameterized by a gyrofactor / and the 
electron gyroradius tl = Ee/{eB). By equating Eq. (2) with Eq. (3), we can estimate 
the gyrofactor as / ~ 1 for both objects, which indicates that the quite efficient 
acceleration has occurred in these PWNe. Interestingly, diffusion of electrons is slow. 
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nearly at the the Bohm hmit of / = 1. The presumed configuration of the toroidal magnetic 
field in PWNe means that the direction of the magnetic field in the post-shock flow should 
be essentially perpendicular to the radial flow direction. Therefore, TeV electrons can be 
conflned quite efficiently in the absence of a radial fleld component. In this regard, the slow 
diffusion is not unexpected for PWNe. If instead diffusion of electrons occurs along a radial 
magnetic fleld line, the magnetic fleld in the post-shock flow must be highly turbulent so 
that the Bohm limit is realized. 

To further investigate diffusion processes, we need to incorporate spatial and temporal 
evolution of the physical parameters of the post-shock flow into morphological and spectral 
modeling. We expect that the next-generation Cherenkov telescopes, Cherenkov Telescope 
Array (CTA), will make it possible to perform spatially- resolved spectroscopy at the TeV 
7-ray band, allowing us to constrai n the spatial and temporal evolution of the PWNe by 



combining X-ray and TeV data (see lVan Etten fc Romanill201ll ). 
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Fig. 12. — Spectral energy distribution for the Suzaku and H.E.S.S. measurements with 
the synchrotron and inverse Compton models. Black thick lines indicate the best-fit results 
with absorbed power-law models. Electron spectral index is assumed to be p=2.3. The 
normalization and the maximum energy of electrons are adjusted to reproduce approximately 
the X-ray data and the TeV data. Black thin lines show the different ambient magnetic fields. 



